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22 Abstract 


23 This study is concerned with the characteristics of storms exhibiting an abrupt temporal 

24 increase in the total lightning flash rate (i.e., lightning jump, LJ). An automated storm 

25 tracking method is used to identify storm “clusters” and total lightning activity from 

26 three different lightning detection systems over Oklahoma, northern Alabama and 

27 Washington, D.C. On average and for different employed thresholds, the clusters that 

28 encompass at least one LJ (LJ 1) last longer, relate to higher Maximum Expected Size of 

29 Hail, Vertical Integrated Liquid and lightning flash rates (area-normalized) than the 

30 clusters that did not exhibit any LJ (LJO). The respective mean values for LJ1 (LJO) 

31 clusters are 80 min (35 min), 14 mm (8 mm), 25 kg m' 2 (18 kg m' 2 ) and 0.05 flash min' 1 

32 km' 2 (0.01 flash min' 1 km' 2 ). Furthermore, the LJ1 clusters are also characterized by 

33 slower decaying autocorrelation functions, a result that implies a less "random" behavior 

34 in the temporal flash rate evolution. In addition, the temporal occurrence of the last LJ 

35 provides an estimate of the time remaining to the storm’s dissipation. Depending of the 

36 LJ strength (i.e., varying thresholds), these values typically range between 20-60 min, 

37 with stronger jumps indicating more time until storm decay. This study’s results support 

38 the hypothesis that the LJ is a proxy for the storm’s kinematic and microphysical state 

39 rather than a coincidental value. 

40 1 . Introduction 

41 The advent of ground-based lightning detection networks in recent decades has 

42 made real-time retrieval of total lightning activity (cloud-to-ground, CG and the intra- 
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43 cloud, IC) available in both high spatial and temporal resolutions. Although there are 

44 uncertainties in the details (Takahashi 1978; Saunders 1993), it is known that 

45 rebounding collisions between graupel and ice crystals in the presence of super-cooled 

46 water is the primary process for thunderstorm electrification (MacGorman and 

47 Morgenstern 1998; Saunders et al. 2006; Emersic and Saunders 2010). Several studies 

48 have documented a temporal co-variability between updraft mass flux, precipitation ice 

49 mass and overall storm depth with the respective total lightning activity (e.g., Goodman 

50 et al. 1988; Carey and Rutledge 2000; Chronis et al. 2007; Deierling and Petersen 2008; 

51 Bruning and MacGorman 2013). Hence, it would be reasonable to suggest that an abrupt 

52 temporal change of the order of a few minutes in the total lightning activity is 

53 considered as a severe weather indicator (“Lightning Jump", LJ, see Schultz et al. 2009; 

54 2011). Studies by Williams et al. (1999), Gatlin and Goodman (2010), Carey et al. 

55 (2009), Schultz et. al. (2009; 2011) and Rudlosky and Fuelberg (2013) document that 

56 statistics such as lead time, probability of detection and false alarm ratio could be 

57 improved based on the use of total lightning as a metric for storm intensity. Nonetheless, 

58 these methods can be hindered by problems related to uncertainties in severe weather 

59 observations at the surface (Trapp et al. 2006; Keene et al. 2008; Schultz et al. 2011). 

60 This study puts forward an original comparison between the convective characteristics 

61 of storms that did or did not exhibit a LJ throughout their lifetime. This evaluation relies 

62 on radar-derived and lightning properties. 

63 2. Data and Methods 
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64 2. 1 Storm Tracking and Clustering 

65 The storm identification and tracking have been performed in real-time utilizing 

66 the Warning Decision Support System Integrated Information tracking system (WDSS- 

67 II, Lakshmanan et al. 2007). A storm “cluster” is automatically identified by the 

68 reflectivity across the -10°C isothermal layer, following a merger of individual WSR- 

69 88D radars. A combination of watershed segmentation and k-means clustering is 

70 employed to identify the storm clusters (Lakshmanan et al. 2009; Kolodziej Hobson et 

71 al. 2012; Cintineo et al. 2014). To complete the storm identification, the algorithm 

72 searches for local reflectivity (Z) maxima where Z > 20 dBZ, then incrementally grows 

73 the area until it is at least 200 km 2 . The storm cluster is then matched with a separately 

74 identified cluster at the next time step (for our analysis, a 1-min time step was used) 

75 using a cost function, where longer-lived cells are given preference in the case of storm 

76 mergers. 

77 Each storm (hereinafter cluster) is described by a geolocated polygon (i.e. 

78 footprint). The cluster’s lifespan is determined as the total time a cluster was identified 

79 and tracked by WDSS-II (Lakshmanan and Smith 2009). The Maximum Vertical 

80 Integrated Liquid (VIL, Greene and Clark 1972) and the Maximum Expected Size of 

81 Hail (MESH, Witt et al. 1998; Cintineo et al. 2012) are retrieved for each cluster for the 

82 duration of its lifetime. Both VIL and MESH have been used as radar-derived intensity 

83 metrics for storm properties such as liquid precipitation, updraft strength and hail growth 

84 (Ambum and Wolf 1996; Witt et al. 1998). As with any proxy, there are caveats that 
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reflect the imperfect representations of severe weather potential and emanate from 
parameters unrelated to the storm dynamics (e.g. distance from the radar, tilted updrafts, 
storm speed etc., Stumpf et al. 2004). To mitigate these effects as much as possible, all 
available radars in the area are used to retrieve these proxies. Five radars over each of 
the three locations are employed, namely, KFDR, KTLX, KVNX, KINX, KSRX for 
Oklahoma, KHTX, KGWX, KBMX, KDHX, KFFC for north Alabama and KLWX, 
KDOX, KAKQ, KCCX, KDIX for DC (radar acronyms from 
https://www.ncdc.noaa.gov/nexradinv/map.jsp). The data for the present study extends 
from 1 April 2013 through 14 August 2013. 

2.2 Total Lightning Activity and the Lightning Jump Algorithm 

This study employs three total lightning detection networks: 1) the Lightning 
Mapping Array (LMA) networks located in central/SW Oklahoma (MacGorman et al. 
2008), North Alabama (Goodman et al. 2005), and Washington D.C (Krehbiel 2008) 2) 
the Earth Networks Total Lightning Network (ENTLN, Liu and Heckman 2010) and 3) 
the National Lightning Detection Network (NLDN, Cummins et al. 1995,2005, 
Cummins and Murphy 2009). 

The LMA networks detect the very high frequency (VHF) radiation emitted 
during the elemental processes that compose a lightning discharge (e.g. the initial 
breakdown, leader propagation and other K-processes, Uman 1987) with a location 
accuracy measured in tens of meters and with a time resolution of 80-100 ps (Thomas et 
al. 2004). The LMA detects both IC and CG flashes although the distinction can be 
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106 


dubious due to limitations in range. The location accuracy is also range-dependent, 


107 however it is relatively constant between -150 km radius from the respective center 

108 (Thomas et al. 2004; Koshak et al. 2004). The following analysis relies on the total 

109 lightning flashes occurring within -120 km of the respective LMA center (Thomas et al. 

110 2003). Lightning flashes are retrieved from the LMAs via grouping at least 10 detected 

1 1 1 VHF radiation sources, using time and space constraints (3 km and 150 ms) between the 

1 12 adjoining points (McCaul et al. 2008). Only flashes that begin within the storm cluster’s 

113 footprint are counted towards the total flash rate. No classification between CG and IC 

1 14 flashes is performed using LMA data. 

115 The ENTLN sensors operate over a wide frequency range, spanning from 1 Hz to 

116 12 MHz. According to Liu and Heckman (2011), electric field waveforms are used in 

117 locating as well as classifying the IC and CG flashes. Multiple strokes (or individual 

118 cloud events) are clustered into a single flash if they are within 700 ms and 10 km of the 

119 first detected stroke. A flash that contains at least one return stroke is classified as a CG 

120 flash, otherwise it is classified as an IC flash. 

121 Since the late 1980s, the National Lightning Detection Network (NLDN, 

122 Cummins et al. 1995; 1998; 2006) has served as the source for many CG lightning- 

123 related studies over the US. The network consists of 113 sensors that combine the 

124 advantages of direction finding and time-of-arrival techniques. The NLDN CG detection 

125 efficiency ranges between 90-95% over the mid latitude continental US, with a median 

126 location error better than 500 m (Cummins and Murphy 2009; Rudlosky and Fuelberg 
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127 2010). Although the NLDN is designed to primarily detect CG flashes, it has been 

128 recently reported that IC flashes are also detected depending on the restrictions applied 

129 to the processed waveforms (peak-to-zero rise time, Murphy and Nag, 2014). 

130 The present study employs the total flash activity (IC+CG) for all lightning 

131 detection systems. Rudlosky and Fuelberg (2013) use a similar methodology for 

132 compiling lightning and radar data. Both NLDN and ENTLN have national (US) 

133 coverage. Nevertheless, for this analysis the respective total lightning activity is 

134 computed only for the clusters that are identified over a radius around where the 

135 optimum LMA operation is ensured. Further detailed comparison (e.g., relative location 

136 accuracy and detection efficiency) between the lightning detection systems lies outside 

137 the scope of this paper. However, their employment is considered as a preliminary 

138 attempt to demonstrate results pertaining to the LJ properties from lightning detection 

139 networks of different technical specifications (e.g., detection efficiency). 

140 The 1-min flash rate is computed by adding all the flashes occurring within the 

141 footprint of the identified cluster. The LJ is objectively identified by Schultz et al. (2009; 

142 2011). This technique uses 14 min of the cluster’s most recent flash rate history. Twelve 

143 of the 14 minutes are considered to calculate the minimum jump threshold that must be 

144 exceeded for a LJ to occur. The remaining two minutes are used to determine whether 

145 the current rate of change in the total flash rate exceeds the LJ threshold. As outlined in 

146 Schultz et al. (2009; 201 1), the algorithm is a 5-step process. These steps are as follows: 

147 1) The total flash rate (f min' 1 ) from the 14 minute period is binned into two minute 
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segments and the total flash rate is averaged (Eqn. 1) 


Fwn = ra(t) 7 (t -« ( i) 


2) The rate of change of the total flash rate (DFRDT, f min' 2 ) is calculated by subtracting 
consecutive bins from each other (Eqn. 2) ^ FR a vg (. 0 — FRav9 ^ ^ _ dfrfT 


( 2 ) 

This results in six DFRDT values (f min' 2 ) 3) The five earliest DFRDT values in time 
are used to calculate the standard deviation (a) of the population 4) If DFRDT > a * a 
and the flash rate is in greater than a given flash rate threshold (FRT) then a FJ has 
occurred. Note that -a- represents a multiplicative factor (i.e. dimensionless) and has no 
relation to the standard deviation (a). In the original studies by Schultz et al. (2009; 
2011) a and FRT were set to 2.0 and 10 f m' 1 . Also note that a is dimensionless (i.e., f 
min' 2 /f min' 2 ). Here we compute the FJ based on a variable a (0.5 - 4, step of 0.5) and 
FRT (5-25 f m' 1 , step of 5 f nr 1 ). The latter is employed in order to define the FJ relative 
strength. For example, a weaker FJ1 would have ofT.O and FRT=10 f m' 1 while a 
stronger FJ1 would be considered as a=2.0 and a FRT=15 f nr 1 5) This process is 
repeated every two minutes as new total lightning flash rates are collected until the 
storm dissipates. 

We note that the above implemented time-window within which the FJ is 
calculated is based on empirical observations of the growth and decay on the convective 
time scale (<10-20 minutes). Had we allowed for longer periods (e.g. 40-60 minutes) 
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168 into the thunderstorm’s lifetime we would likely have missed the occurrence of the first 

169 LJ and potentially severe weather occurrence. This is why we’ve empirically tested this 

170 algorithm with over 700 storms in multiple storm environments to help understand the 

171 variability of the algorithm (Schultz et al. 2011). The choice of the 2 *g (i.e. a=2) in 

172 Schultz et al. (201 1) is simply a benchmark to which this study is not tied to. 

173 3. Analysis and Discussion 

174 3.1 Data and Quality control 

175 As WDSS-II tracks clusters independently of the respective total lightning 

176 activity, the number of the identified LJ0 and LJ1 clusters is considerably different. For 

177 instance more than 2,000 clusters are classified as LJ0 at a=2.0 f m' 2 and FRT=10 f m' 1 

178 whereas less than 200 are classified as LJ1 for the same a and FRT values. To ensure a 

179 comparable sample size and improve the representativeness of the data, we report on the 

180 LJ0 and LJ1 clusters that exhibit sustained total lightning activity for more than 95% 

181 during their lifespan (e.g., if a cluster is tracked for 100 minutes, the cluster must exhibit 

182 total lightning activity greater than zero for at least 95 minutes). This quality constraint 

183 (QC1) may classify a slightly different number of clusters depending on the employed 

184 lightning detection system. An additional quality constraint (QC2) is applied to the 

185 clusters that start or end at a flash rate that is notably higher than zero (set to >10 f m' 1 ). 

186 Typically, these cases represent merging or splitting clusters or clusters that 

187 entered/exited the effective radius of the LMA with high flash rates. QC2 also takes 

188 care of potential problems with MESH/VIL repetitiveness due to distance from the 
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189 radar. Given that the study explores aspects such as the storm duration, the clusters that 

190 failed to conform to QC2 are omitted from the analysis. Figure 1 illustrates examples 

191 from a tracked cluster that exhibits a problematic tracking (e.g. cluster entering the area 

192 with already high flash rates, Fig. la), a normal tracking (i.e. comply with both QC1 and 

193 QC2, Fig. lb), a LJO cluster (Fig. lc) and a LJ1 cluster (Fig. Id) of comparable flash 

194 rates. 

195 All three lightning detection systems indicate that the number of LJ1 clusters 

196 decreases as the a and FRT values increase (i.e. fewer clusters at higher a and FRT 

197 values, Fig. 2b, d, and f). Unlike the LJ1, the number of tracked LJO clusters increases as 

198 the values of a and FRT increase (Fig. 2a, c, and e). The latter should be expected since 

199 a LJO at e.g. a=2.0 and FRT=10 f m' 1 will also not exhibit LJ at higher a or FRT values. 

200 3.2 The Autocorrelation function of LJO and LJ1 Flash Time Series. 

201 Autocorrelation is an essential tool for describing the independence of sequential 

202 values in a time series. A slow (fast) decaying autocorrelation function with time (i.e. 

203 lag) indicates a consistent (random) behavior of the variable under consideration 

204 (Bowerman and O'Connell 1979). For example, a slow-decaying autocorrelation 

205 function of lightning activity time-series would signal a coherent behavior in the storm’s 

206 updraft speed and volume (e.g., Schultz et. al. 2009; 2011; Schultz et al. 2014). 

207 Consequently, autocorrelation can elaborate on whether the presence of a LJ relates to a 

208 numerically random increase in the total lightning activity or points to a more persistent 

209 feature of the storm's dynamical evolution. The autocorrelation function is computed for 
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210 the flash rates of LJO and LJ1, by introducing a time lag that ranges from 1 to +N/2 

211 minutes, where N is the number of 1 -minute intervals during which the cluster is tracked 

212 (i.e., lifespan). The lag at which the Pearson correlation is reduced below the 95% 

213 significance level denotes the “ e-folding ” time.. Figure 3 illustrates the average e-folding 

214 times for the LJO and LJ1 clusters for different a and FRT values. The corresponding 

215 results (Fig. 3) show longer e-folding times for the LJ1 clusters. For example, the e- 

216 folding times for the LJ 1 at a=2.0 and FRT=15 f m' 1 are computed as ~12 min for LMA, 

217 12.7 min for the ENTLN, and 1 1.5 min for the NLDN. Conversely, the e-folding times 

218 for the LJO for the same a and FRT values are consistently less than -4.0 min for all 

219 three lightning detection systems and any given a and FRT value. Moreover, the fact 

220 that the e-folding times for LJ1 clusters increase as both a and FRT values also increase, 

22 1 illustrates a key observation that emphasizes the non-redundant numerical role of both 

222 variables a and FRT in the LJ algorithmic implementation (Schultz et al. 2009; 2011). 

223 3.3 Comparison of storm severity potential and physical characteristics between 

224 the LJO and LJ1 clusters. 

225 The previous section studied the LJO and LJ1 clusters exclusively from the 

226 standpoint of the flash rate temporal variation. This section explores the mean values of 

227 storm attributes derived from WDSS-II. As Fig. 4 demonstrates, the LJ1 clusters exhibit 

228 a longer lifespan than the respective LJO, and this observation is consistent throughout 

229 the three lightning detection systems and all a and FRT values. For example, for a=2.0 

230 and FRT=15 f m" 1 , the average lifespan is 80 min, whereas the respective LJO lifespan is 
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231 approximately 35 min. Similar behavior is evident for the mean flash rate (normalized 

232 by the cluster’s footprint area, f m' 1 km' 2 , Fig. 5), MESH (Fig. 6) and VIL (Fig. 7) 

233 values. 

234 In particular, Fig. 5 indicates that on average, the LJ1 clusters exhibit -4-5 times 

235 higher flash rates than the respective LJO. For instance, the average LJ1 flash rates for a 

236 =2.0 and FRT=15 f m' 1 are -0.054 f m' 1 km' 2 as opposed to -0.015 f m' 1 km' 2 for the 

237 LJO, an observation that is also consistent across all networks. In turn, the MESH values 

238 for the LJ1 clusters range from -1 1-18 mm whereas the LJO corresponding values range 

239 from -6.5-10 mm (Fig. 6). Likewise, the mean values of VIL are - 18 kg m' 2 for the 

240 LJO and -25 kg m' 2 for the respective LJ1 (Fig. 7). As also highlighted in Section 3.2, 

241 higher flash rates, larger MESH and VIL values (Figs. 5-7) are found with increasing a 

242 and FRT thresholds. One could argue that it would be expected to have higher 

243 magnitudes of weather severity proxies (e.g., MESH, VIL etc.) with higher flash rates. 

244 Nevertheless, the previous results also suggest that it is not only the flash rate (i.e., FRT) 

245 that exhibits a fundamental physical tie to storm intensity but also its temporal evolution 

246 (i.e., a). The above results are also in agreement with the findings by Rudlosky and 

247 Fuelberg (2013). 

248 3.4 LJ strength and storm decay time 

249 The results shown in Fig. 4 support that the LJ1 clusters with larger a and FRT 

250 relate to storms with longer durations (Fig. 4). Approaching this from a different 

251 perspective one can raise the following question: “Does the strength of the final LJ 
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252 occurrence relate to the remaining lifespan of the cluster?” To address this question the 

253 time elapsed from the last LJ occurrence to the last time-step that a cluster is tracked is 

254 computed in minutes. Arguably, the results in Fig. 8 corroborate that both a and FRT 

255 play a role in the storms’ remaining duration which shows to increase from around 30- 

256 35 min for LJ1 of a = 1.0 and FRT=10 f m' 1 to over 45-60 min for higher a and/or FRT 

257 values. 

258 4. Conclusions 

259 The observations herein indicate that the presence of LJ has implications for the 

260 storm dynamics, intensity and evolution. The e-folding times are lower for the LJ1 

261 clusters. For example the e-folding times for the LJ1 at a =2.0 and FRT=15 f m' 1 are 

262 computed as -12 min for LMA, 12.7 min for the ENTLN, and 1 1.5 min for the NLDN. 

263 Conversely, the e-folding times for the LJ0 for the same a and FRT values are 

264 consistently less than -4.0 min for all three lightning detection systems. Through the 

265 enhanced updraft hypothesis, these findings indicate that the presence of a LJ signals the 

266 storm’s ability to sustain convection. 

267 The study also documents that LJ1 clusters last longer and exhibit higher flash 

268 rates (area-normalized), MESH and VIL values, further corroborating previous studies 

269 that also suggest that the temporal total lightning variability is a dependable proxy for 

270 severe weather risk assessment (Williams 2001; Schultz et al. 2009; 2011; Rudlosky and 

271 Fuelberg 2013). In addition, the MESH values for the LJ1 clusters range from -11-18 

272 mm whereas the LJ0 respective values range from -6.5-10 mm (Fig. 5). The mean 
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273 values of VIL are ~ 18 kg m' 2 for the LJO and ~25 kg m' 2 for the LJ1 clusters. 

274 The results throughout this analysis consistently suggest that there is no 

275 redundancy in the role of a and FRT in the LJ numerical implementation. This is shown 

276 by the increasing magnitudes of the implicated variables (e.g. e-folding time, MESH, 

277 flash rate etc., see Fig.2-7) for LJ1 clusters increase as both a and FRT values also 

278 increase. Finally, the study offers further evidence that the presence and temporal 

279 coincidence of a LJ could be viewed as a proxy of the storm’s expected dissipation. 

280 Typically, these values range between 20-60 min depending on the LJ strength with 

28 1 stronger jumps indicating more time until storm decay. 

282 Ongoing efforts explore the value of the LJ as a component in the operational 

283 severe weather watch/warnings issuance (Schultz et al. 2014). The upcoming Geo- 

284 stationary Lightning Mapper (GLM) onboard the GOES-R mission (Goodman et al. 

285 2013) will provide continuous monitoring of total lightning activity across the Western 

286 Hemisphere. GLM will provide even greater detail on the linkage between temporal 

287 lightning variability and the storm evolution over areas where currently related 

288 information, including radar, is limited or nonexistent. Importantly, GLM will provide 

289 continuous coverage of total lightning over a large domain to evaluate this study on the 

290 global scale. 
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495 Figure Captions 

496 Figure 1: Tracked cluster that exhibits a problematic tracking (e.g. cluster entering the area with already high flash rates, 

497 Fig. la), a normal tracking (i.e. comply with both QC1 and QC2, Fig. lb), a LJ0 cluster (Fig.lc) and a LJ1 cluster (Fig. Id) of 

498 comparable flash rates. 

499 

500 Figure 2: The identified number of LJ0/LJ1 clusters as a function of FRT (x-axis, f nr 1 ) and (X (y-axis,), LMA- 

501 a/b, ENTLN-c/d and NLDN-e/f 

502 Figure 3: Mean e-folding time (min) for LJ0/LJ1, as a function of FRT (x-axis, f nr 1 ) and (X (y-axis,), LMA- 

503 a/b, ENTLN-c/d and NLDN-e/f 

504 

505 Figure 4: Mean life-span (min) for LJ0/LJ1, as a function of FRT (x-axis, f nr 1 ) and (X (y-axis,), LMA-a/b, 

506 ENTLN-c/d and NLDN-e/f 
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Figure 5: Mean area-normalized flash rate (f m" 1 km' 2 ) for LJ0/LJ1, as a function of FRT (x-axis, f m" 1 ) and CL 
(y-axis,), LMA-a/b, ENTLN-c/d and NLDN-e/f 

Figure 6: Mean MESF1 (mm) for LJ0/LJ1, as a function of FRT (x-axis, f m" 1 ) and CL (y-axis,), LMA-a/b, 
ENTLN-c/d and NLDN-e/f 

Figure 7: Mean VIL (kg m" 2 ), for LJ0/LJ1 as a function of FRT (x-axis, f m" 1 ) and CL (y-axis,), LMA-a/b, 
ENTLN-c/d and NLDN-e/f 

Figure 8: Time elapsed until the storm dissipation for LJ1 (min) (LMA-a, ENTLN-b and NLDN-c) as a 
function of FRT (x-axis, fm' 1 ) and CL (y-axis,). 
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543 Figure 2: The identified number of LJ0/LJ1 clusters as a function of FRT (x-axis, f m" 1 ) and CL (y-axis,), LMA- 


544 a/b, ENTLN-c/d and NLDN-e/f 

545 

546 

547 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 

560 

561 

562 


30 of 45 















564 Figure 3: Mean e-folding time (min) for LJ0/LJ1, as a function of FRT (x-axis, f m" 1 ) and CL (y-axis,), LMA- 
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586 Figure 4: Mean life-span (min) for LJ0/LJ1, as a function of FRT (x-axis, f m" 1 ) and (X (y-axis,), LMA-a/b, 


587 ENTLN-c/d and NLDN-e/f 
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607 Figure 5: Mean area-normalized flash rate (f m" 1 km' 2 ) for LJ0/LJ1, as a function of FRT (x-axis, f m" 1 ) and (X 
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630 Figure 6: Mean MESH (mm) for LJ0/LJ1, as a function of FRT (x-axis, f m" 1 ) and (X (y-axis,), LMA-a/b, 
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651 Figure 7: Mean VIL (kg m" 2 ), for LJ0/LJ1 as a function of FRT (x-axis, f m" 1 ) and (X (y-axis,), LMA-a/b, 
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673 Figure 8: Time elapsed until the storm dissipation for LJ1 (min) (LMA-a, ENTLN-b and NLDN-c) as a 


674 function of FRT (x-axis, f m' 1 ) and (X (y-axis,). 
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